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CHAPTER 1 . INTRODUCTION 


Pressure vessels which are made of filamentary com- 
posites have had a profound impact in applications where 
weight efficiency is a factor. In particular, filamentary 
composite pressure vessels with a spherical geometry are 
used in applications of space shuttle tankage [1] and for 
storage of various fluids [2]. 

Spherical storage pressure vessels usually consist of 
a thin metallic bladder onto which is wound high-strength 
fibers bound by an organic matrix. The bladder, which acts 
as a liner to prevent leakage, is usually fitted with one 
or more fill tubes for introducing the fluid to be pressur- 
ized. As shown in Figure 1, the spherical shape allows 
filament winding patterns which result in a laminate which 
may be considered as pseudo-isotropic. That is, the 
pseudo-isotropic (transversely isotropic) condition is met 
if the angles between lamina in a laminate satisfy the 
relation [3]: 

TT 

Y • (1) 

M 

where M - 3,4,3... 

This wrap angle ensures that the stiffness matrix will be 
invariant for the 0-41 plane shown in Figure 2 and the 
composite may be modeled as a transversely isotropic 
material . 
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FIGURE 1, Cui^way of typical filamentary composite 
spherical pressure vessel [7] 
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A number of closed-form solutions have been developed 
[3,4] which predict the stress-strain behavior of a compos- 
ite shell assuming a pseudo-isotropic material and an 
elastic-plastic bladder response. These solutions predict 
short-term failure pressures of the vessel, and correlate 
well with experimental data [4,5,6]. Finite -element 
techniques are also used [2,7] to model discontinuities 
near the fill tubes. 

Although recognized as a significant problem, little 
is known about the mechanisms which lead to the loss of 
long-term structural integrity of filamentary composite 
spherical pressure vessels. Some possible causes of this 
time -dependent loss of structural integrity are environmental 
degradation, matrix crazing, and manufacturing processes. 

The static fatigue susceptability in composites is 
well-known, but it has not been treated successfully from 
the standpoint of mechanism identification. Therefore, in 
the case of long-term, high-pressurization storage, high 
factors of safety are used to condensate for uncertainties 
regarding creep rupture under conditions of sustained 
pressurization. The resulting low design allowables then 
detract from the weight-saving attributes of composites. 

This problem has heretofore been dealt with by a statistical 
analysis of manufacturing variations and experimental data 
[ 8 ]. 
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The study of the time -dependent behavior of filamentary 
composite pressure vessels is further complicated by the 
lack of experimental data. Some data t?ere collected at 
Lawrence Livermore National Laboratory on S-glass/epoxy 
composites [9], but most of the test specimens were destroyed 
and the test data invalidated due to a major earthquake in 
Northern California in 1980. 

This study examines the time -dependent response of 
filamentary composite spherical pressure vessels by 
considering the linearly vis^:oelastic character of the 
matrix in the filamentary composite laminate. This 
characterization is important because the composite material 
over a period of time becomes more anisotropic, thus in- 
creasing the failure potential of the system. Therefore, 
this material model which includes the time -dependent 
properties of the matrix in a composite laminate is used 
to predict burst pressure and strain under conditions of 
sustained internal pressurization of a spherical pressure 
vessel. The lamina and laminate properties in the analysis 
are initially determined using micromechanical and macro- 
mechanical techniques. Viscoelastic material properties for 
the matrix are experimentally determined and included in the 
constitutive equations for the composite. These time- 
dependent expressions are transformed into the Laplace domain 
and solved using the maximum strain theory of failure. The 
solution is then inverted back into the time domain using the 
method of collocation, and burst pressure and critical strain 
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of the composite spherical pressure vessels are expressed 
as a function of time. Finally, an application is given 
for a Kevlar/epoxy system, where the solution is numerically 
evaluated. The burst pressures for this system are predicted 
as a function of time, and critical strain is computed for 
a constant internal pressurization. The effects of 
constituent volixne fractions and shell thickness on failure 
potential are also investigated. Thus, an analytical method 
is outlined and a numerical application is given to predict 
the time -dependent response of a filamentary composite 
spherical pressure vessel. 


CHAPTER 2. 


MATERIAL MODEL 


Fllai&entary composite spherical pressure vessels general- 
ly have filaments which are wound into a la'^ilnate considered 
as pseudo-lsotroplc (or transversely isotropic) . In this 
analysis, the filamentary composite pressure vessel is as- 
sumed to be transversely isotropic in the 6 - 0 plane, as 
shown in Figure 2. 

The spherical geometry of pressure vessels unde,r consid- 
eration is uniform, except in the areas surrounding the fill 
tubes. However, the model used in this analysis assumes a 
uniformly spherical geometry and does not include the dis- 
continuities in the fill tube regions. 

It has been shown that the bladder which lines the inner 
surface of the composite shell has a minimal effect in the 
burst pressure prediction of the pressure vessel [4, 10]. 
Therefore, the structural model used in this analysis does 
not include the effect of the bladder. Models have been 
developed which consider the elastic-plastic behavior of 
the bladder in failure prediction [3]. 

This model will be used in the following viscoelastic 
analysis to predict burst pressure and to calculate critical 
strain. It is further elaborated upon in the application 
to a Kevlar/epoxy spherical pressure vessel. 
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CHAPTER 3. EQUATIONS OF ELASTICITY 

For a spherical geometry with uniformly symmetrical 
loading, all of the shear stress and strain components 
vanish. For the special case of a transversely isotropic 
material system, the general governing equations of elastic- 
ity may be further simplified. The following discussion 
derives the equilibrium equations and kinematic equations 
in spherical coordinates for a sphere transversely isotropic 
in the r plane. 


3.1 EQUILIBRIUM EQUATIONS 

Neglecting body forces, the following general equilib- 
rium equations relative to a spherical coordinate system are 
[ 111 : 


aa 


^ + I 


aa 


ao 


M + 3^ , + _£214 ® ^2a) 


ar 


r a<f rsin<j) ae r ^ rsin()> 


a a 


da 


re ^ 1 


da 


ee , 3 , 2cos(j) ^ _ n 

+ ^ CJa. - 0 


R Rfh 

3r r 30 rsin0 36 r rsin0 ^ 

Bo, 


(2b) 


Bo , 3a 1 

EL + ± 1 

3r 


r0 , cos(t> 

* 0 . j » 


r 30 rsin0 36 rsin0 


; - “eel ' ° 


(2c) 


where are stress components in the 6,0,r coordinate sys- 
tem. For uniform loading conditions, all of the shear stress 


components vanish, i.e., a^^ = a^^ = aj,g 
Therefore, Equations (2) reduce to: 


0r 


06 ^60 


0 
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rsin4i 30 


3o 1 

-^ + i (2o - 0 

3r r " ' 


- ‘=ee’ = ° 


(3a) 


(3b) 


(3c) 


For a transversely isotropic material system, the tangential 
stresses, and are equal. Equations (3) thus can be 
written: 



r d4> 


da 

rsinG d9 


09 


= 0 


2 

+ - (0„ - o.„) - 0 

dr r “ 


(4a) 


(4b) 


(4c) 


Equation (4c) is then the only equilibrium equation of 
interest . 


3.2 KINEMATIC EQUATIONS 

For linear elasticity, the kinematic equations in spher- 
ical coordinates for any geometry are given by [ 11] : 


•(j)(|) 


r 3(|) r 


(5a) 
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r- « 1 i 

rsin4> 3^ 


rr 


3r 


'64) 


i ^ cot4> + 1 — a ) 


rsln4> 36 


r 34) 


■<^r 


i ( i 4 ^ ) 


r 3(^ 


3r 


^ . 1 r !lis+ He ^ 

'r 6 ^ ~ 

2 3r rsin4> 36 r 


(5b) 

(5c) 

(5d) 

(5e) 

(5f) 


vrtiere are the strain components and are the displace- 
ments in the 6,(p,r coordinate system. 

For a condition of uniform loading, all shear strain 
components vanish and Equations (sd) through (5f) may then 
be eliminated. Also, displacement components u^ amd Ug are 
zero and Equations (5a) through (5c) reduce to: 

" (6a) 


■60 


rr 


r 

lir 

r 

dr 


(6b) 

(6c) 
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For a material transversely isotropic in the r plane, 
^66 “ a*'*! Therefore, the final form of the 

kinematic equations is: 


^ee 


u 

r 


e 


rr 


du 

dr 
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CHAPTER 4. MACROMECHANtCAL BEHAVIOR OF A LAMINA 


A lamina In a filamentary composite is an assemblage 
of fibers embedded in a matrix. It is necessary to describe 
the macromechanical behavior of a lamina, l.e., the behavior 
when averaged properties are considered, in order to under- 
stand the laminated structure. In this discussion the 
coordinate system 6,(|),r is aligned with the principal material 
directions 1,2,3 of the lamina. 

The generalized Hooke's law relating stresses to strains 
can be written as: 

'^ij “ ^ijkl^kl i.j.k,l » 1,2,.. .,6 (8 ) 

where are stress components, is the elasticity ma- 

trix, and are the strain components. Although the elas- 
ticity matrix has 36 terms, it is symmetric and only 21 of 
the terms are independent. 

The 21 constants in the elasticity matrix are necessary 
in order to characterize a general anisotropic material be- 
cause there are no planes of symmetry. The case of a lamina 
which is transversely isotropic in the 3 plane yields prop- 
erties which are non -directional in the 1-2 plane. The 
stress-strain relationship is based upon only five independ- 
ent constants and is given by [12 ] : 


12 
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(9) 


Letting subscript 1 “ 0» subscript 2 • <p, and subscript 3 ■ r 
in spherical coordinates, the stress-strain relations become: 


^^00 


<=00 <=04> <=0r 


^00 

^<P<P 
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CD 

-o- 

o 

CD 

CD 

o 

CD 



^rr 


|^0r <=0r <=rrj 


Jrr^ 


Again, due to loading symmetry, 
constitutive equations reduce to: 

^00 " ^^60 ^00 ^0r ^rr 

and 


( 10 ) 



Therefore, the 


(11a) 


^ ^0r ^00 ^rr ^rr 


(11b) 


CHAPTER 5. MICRC»*ECHAHICAL BEHAVIOR OF A UMIHA 


Micromechanics is the study of coi^osite behavior «7here 
the Interaction of the constituent materials Is examined In 
detail. This study Is based on a lamina consisting of two 
constituents, the fiber and the matrix. In the coordinate 
system, L,T,r , the longitudinal direction, L, Is along 
the length of the fiber, and the transverse direction, T, 

Is across the width of the fiber and Is orthogonal to L. 

The radial direction, r. Is perpendicular to the L*T plane. 

The constitutive relationship for an Individual lamina 
In the 6-0 plane may be expressed by the matrix given In 
Equation (10). If the coordinate axes (L,T,r) are oriented 
In the principal material directions and assumed to be 
transversely Isotropic In the r-plane perpendicular to the 
fiber direction, the elasticity matrix is [3 ]: 


[C] 


^LL 

^LT 


0 

0 

0 

^LT 

^TT 

^Tr 
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S.r 
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Crr 

0 
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®Tr 
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0 
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0 
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0 
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( 12 ) 
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Because of isotropy in the T-r plane: 


^Lr “ ^LT 


rr 




°Lr - 


and 


G . ~ 

Tr 


The corresponding lamina compliance matrix is: 



Inversion of the lamina con^liance matrix yields the 
elasticity matrix. Equation (12), and the following 
relationships (13J: 
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(13a) 

(13b) 

(13c) 

(13d) 


(14a) 
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C-rr “ ^1 “ ^LT 

^LT " '’TL 

^Tr " ^^TT '*' ^LT '^TL^ 

Gxr - 

wh^»re 

A ■ [ (1 + v^ip) (1 ” ^TT ” ^^LT '^TL^ ^ ^ 
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( 14b) 


The lamina engineering constants (El* E^, Vj^^* '^xl* '^TT' 
expressed in terms of the elastic properties 
(E, V, G) and voliMne fractions of the fiber and matrix 
of the composite material. Micromechanical techniques 
are used which are based upon the assumption that the 
properties of the fiber and matrix ar« isotropic and are 
equal in tension and compression. 

5.1 DETERMINATION OF E^ 

For a volume of material loaded in the longitudinal 
direction, basic mechanics of materials says that 




A L 
L 


(15) 


where L is the specimen length and is the longitudinal 
extensional strain. 
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The stresses in the fiber and the matrix are given by 




e 


m 



(16a) 

(16b) 


where Ef is the fiber modulus and E is the matrix modulus, 
r m 

The average stress a acts on a cross-sectional area A, 

acts on the cross-sectional area of the fiber, A£, and Qjjj 

acts on the cross-sectional area of the matrix, A^. There- 

m 

fore, the resultant force on the composite material is: 


P = OA - aj Af + A^ 

knowing that 

a = 

substituting Equa‘-ions (16) and (18) into Equation (17) 
and rearranging yields : 




+ E 


m 



(19) 


Defining the volume fractions of the fiber and matrix, 
respectively, 

''f - ~r 

and substituting Equations (20) into Equation (19) results 
in the expression: 


E 


L 


Ej Vj 


+ E. 


m 


V 

m 


( 21 ) 
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Equation (21) is the "rule of mixtures" expression for the 
elastic modulus in the direction of the fibers. 


5.2 DETERMINATION OF E.j. 

If a stress is applied in a direction transverse to 
the fibers, the strains in the fiber and in the matrix are 
given by: 


e 


f " 



(22a) 


e 


Ot 


Let W be the transverse dimension of the specimen 

fore, on an average basis, acts on v^W and 

V W. Thus, the total transverse deformation is 
m 


(22b) 

. There - 
acts on 


or 


v^W 


+ V W 
t m 


E 

m 


(23) 


where is the transverse extensional strain. 
Substituting Equations (22) into Equation (24) gives 


(24) 


'f E. 


+ v_ 


m 


(25) 


m 


Knowing that cr-p ® E^e.p , (26) 

and subr.ituting this relationship into Equation (25) yields 
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Om 

T 


% (Vf 




(27) 


Solving for E,j, from Equation (27) provides the expression 
for the transverse modulus: 


E^E 
f m 


v_E^ + VrE 
mi I m 


(28) 


5.3 DETERMINATION OF AND 

The major Poisson's ratio, can be determined by a 

method similar to that used to determine E^. For a volxune 
of material loaded in the longitudinal direction, the 
Poisson's ratio is defined as: 



For a deformation in the transverse direction: 


AW = — W e,j. 


(30a) 


or 


AW W 


(30b) 


Tbe total deformation is also given by 

AW = AW^ + AWjj^ ( 31 ) 

where W^ is the fiber width and W_ is the matrix width, 
r m 

On the average basis, acts on v^W and e acts on v W. 

I I m m 
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Therefore, 


and 


* Wv£e£ ® 


AW^ • Wv„e„ ® 

m mm m m L 


(32a) 


(32b) 


where Vf is the Poisson's ratio of the fiber and is the 
£ m 

Poisson’s ratio of the matrix. 

Then , 

AW = W + VfVfCi,) (33) 

Substituting Equation (30) into Equation (33) gives 


''lt + '"£''£^ 

which reduces to the expression: 


LT mm ft 


(34) 


(35) 


Equation (35) is the "rule of mixtures" expression for the 
major Poisson's ratio. 

Because of S 3 rmmetry requirements of the lamina compliance 
matrix, can be obtained from the relationship: 


'TL 


XT 


( 36 ) 


5.4 DETERMINATION OF v.p.j, AND G^.j. 

Although there is not general agreement about the 
prediction of transverse lamina properties, reasonable 
predictions may be expressed as [14,15]: 
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1 + 


V, 


TT 


+ v_v_ 
£ r mm 


( 


1 - V. 


m 


'^m^LT ( gS) 


(37) 


G = ^ (4-it) t tlN ^ 4N ^ 

^^2 4 (4-tt) N + ir 

where 

Gf(TT+4v^ + G^(7T-4vp 
G^(Tr-4v^) + G^(tt+4v^) 


(38) 


( 39 ) 


CHAPTER 6. PROPERTIES OF THE LAMINATE 
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Laminate properties may be determined by averaging 
properties of the Individual lamina [12] • The laminate 
stiffnesses C^g, Cgg, Cg^) are obtained by approximat- 

ing the various lamina orientations In the spherical geometry 
as a continuous function of ot, shown In Figure 3. An 
average stiffness Is then computed by Integrating the lamina 
stiffnesses which have been transformed from the (L,T,r) 
coordinates to the (l,2,r) coordinate system, also shown In 
Figure 3. These averaged elastic properties are derived 
with the assumption of constant strain In the laminate, 
Independent of the orientation. 

The average laminate stiffness In the coordinate 

system Is calculated from: 

‘^ee ■ ? (4oa) 

O 

O 

Sr (40c) 

o 

and Cj.j. = (40d) 

These equations may be used to determine the laminate elastic 
stiffness for any thickness. If the Individual lamina Is 
thin compared to its width. 
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Using the method outlined by Gerstle [ 3 J, the 
following equatione are obtained: 

C,, - , i ^ 

S* - + -P + i (3C^^ . 2G^j) 

C__ = — fr. 4- r* \ 


er “ 2 


Assuming that the transversely isotropic 
in the 0 - (ji plane 


” ^06 


G = ^90 ~ ^ 9 ( 

94> ^ 


^4>r ■ ^0r 


condition is met 


By substituting Equations (14 ) into Equations (41 ) . the 
laminate stiffnesses may be obtained in terms of the 
elastic constants of an individual lamina: 


*'6* . ( ^^■''tt^)^t+C1-Vt^v«. )E_+4>i_ nj... 1 .®LT 


TT ^ ^TL ^ ^TL ^ ^ ^ ’ 


(43a) 


(43b) 


‘^er* 3 ‘^lt + 


(43c) 


Ej. A(1 -Uj_j,U3,^) 


(43d) 
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then obtained from: 


ee 


2 Cer 3 


2 C 


00 rr 
2 


- C 2 

er 


0r 


(44) 


“rr 

rr 

^00 ^04) 

(45) 

\) s 

^04 

^rr ’ ^0r 

(46) 

6(t) 

o 

CD 1 

CD 1 

.... . 
"rr - "er 


O 

CD 

^^00 ■ ^0*3 

(47) 

^dr “ 

C o rt 

C „ - C._/ 


Finally, the degree of anisotropy of the laminate is 
obtained from the relationship: 

3 . % 

2 


" = I (^) 


(48) 


rr 


It can be noted that if n - 1.5 the composite material is 
isotropic. The magnitude of n is affected by factors such 
as wrap angle, fiber and matrix content, and, as will be 
demonstrated herein, the time -dependent properties of 
the matrix material . 


CHAPTER 7. 14AXIMUM STRAIN THEORY OF FAILURE 


Prediction of the burst (failure) pressure of a filamen- 
tary composite spnerlcal pressure vessel with a metallic 
bladder Is a function of the elastic-plastic behavior of 
the metallic bladder as well as of the composite material, but 
it has been demonstrated [4,10] that the behavior of the blad- 
der has relatively little effect on the failure prediction of 
the composite. Therefore, the effect of the bladder on the 
burst pressure or circumferential strain of the composite 
sphere will not be considered In the following formulation. 

Prediction of the burst pressure requires a determina- 
tion of the maximum values of stress and strain at the Inner 
surface of the sphere . The maximum strain theory of failure Is 
chosen for this study for the following reasons: 

1) Considerable experimental data has shown the general 
applicability of the maximum strain criterion for predicting 
vessel failure. 

2) The published values for composite maximum strain 
have been shown to be essentially Independent of vessel 
configuration or fiber volume fraction. 

3) The calculated maximum strain includes the contribu- 
tion of radial stress. This is Important in predicting 

the behavior of thick -walled vessels. 
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7.1 


BOUNDARY CONDITIONS 
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Boundary conditions for the composite material portion 
of the pressurized structure are: 

<3.j. - “P at the inner surface, r ■ ri (49a) 

*^rr " ® outer surface, (49b) 

where o^j. is the radial stress component and p is the inter 
nal pressure. 

These boundary conditions and the symmetry of the 
spherical geometry greatly simplify the analysis and pre- 
clude the necessity of a three-dimensional analysis. 


7 . 2 FAILURE PREDICTION 

Within the composite, Equations (4), (7), and (11) 

yield: 


d^u . 2 Crr du 


’ r- 4- 

" d? r dr 


r ^ 2 (c,e - Cee - ^ * <> 


(50) 


The solution to t;iis second order differential equation is: 


u » ar*'”^ + br 


,-n-% 


(51) 


where a and b are constants given by: 



and f and g are constants given by: 
f - (n - %) Cj.j.-H 2Cg^. 

g - (n + %) Cj.j.- 2CQr 


(52) 

(53) 

(54) 

(55) 
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The Tnaxirauin strain prediction is found from : 
a r["*7^+ b 


( 56 ) 


where is the ultimate composite tensile strain. It is 
reasonable to assume that the maximtun strain in the structure 
will be no greater than the ultimate tensile strain in the 
fibers. Experimental results from biaxially loaded pseudo- 
isotropic cylinders [ 3 ] have shown failure strains of 

about 0.75 cr, where ■ — if the fibers remain linear 

r t E£ 

to failure. Then, 


•'e* O-” E7 


(57) 


where ■ Fiber modulus of elasticity 
F£ - Fiber fracture strength 


The maximum strain theory of failure predicts that vessel 
failure occurs when the cirexunferential strain in the compos- 
ite is equal to or greater than the composite ultimate 
tensile strain, K^. This failure criterion may be expressed 
as : 


e 


00 




(58) 


The cosmos ite strain, e.., is found by evaluating 

0 0 

kinematic Equation (5b) for the displacement given in 
Equation (6b). 

Thus , 


(-)+ 6^" 


g(B^"-l) f 
where 8 is the ratio r^ / r^. 


( 59 ) 
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Letting the pressure, p, equal the burst pressure, p£, and 
comparing and the following equation can be used 
to determine when composite failure will occur: 


Ke 1 




(-)+ 

f 


If failure occurs when Kg equals tgg, 
burst pressure Is given by: 


pf ■ Kgg 


(6^^-l) 


g 

£ + 6 


2n 


(60) 

the expression for 

(61) 


With the maximum strain theory of failure, Equations 
(60) end (61) are used to predict composite failure for a 
trt^nsversely Isotropic spherical pressure vessel . These 
equations form the basis for the following derivations which 
consider the time -dependent behavior of the composite. 




CHAPTER 8. DETERMINATION OF VISCOELASTIC MATERIAL PROPERTIES 


A stress analysis of viscoelastic materials requires the 
use of mechanical properties which characterize the time- 
dependent behavior of the materials. These material proper- 
ties must be experimentally determined. In this discussion of 
the properties of a linearly viscoelastic material, only 
the epoxy matrix of the composite is considered to exhibit 
t ime -dependent behavior . 

The determination of the relaxation im>dulu8, E^(t), is 
fundamental to the characterization of a linearly viscoelastic 
material. In a standard creep test, a constant uniaxial 
stress 0 ^ is applied and the time -dependent uniaxial strain 
is measured. The time -dependent strain is proportional to 
and may be written as 

^(t) - (62) 

where is Che creep coaq>Iiance of the material. The 

creep coo^llance is the strain per unit of applied stresti, 

^d describes uniquely the stress-Ptrain behavior of a linear- 
ly viscoelastic material. 

8.1 LAPLACE TRANSFORMATIONS 

If Oq is set equal to unity and the Laplace transform 
of Equation (62) is taken, the following result is obtained: 

e(s) “ ^_(s) (®3) 

m 
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%fhere s is the Laplace parameter and the Laplace transform 
is defined by 
00 

f(s) = f (64) 

o 

If a standard stress relaxation test is now considered, 
wherein a uniaxial displacement is applied resulting in a 
normal strain which is held constant, the corresponding 
time-* dependent uniaxial stress is related to the constant 
uniaxial strain through the relationship 

cr(t) = (65) 

If is set equal to unity and the Laplace transform of 
Equation (65) is taken, the following result is obtained; 

a(s) = Ejjj(s) (66) 

In the Laplace domain the "associated elastic" stress-strain 
relationship may be written 

a(s) = Ejjj(s) e(s) (67) 

and the corresponding strain-stress relationship may be 


written 
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where e(s) and J(s) are the associated elastic modulus 
and compliance, respectively. If conditions of the standard 
uniaxial creep test, i.e., 

E(t) = Jjjj(t) (69) 

where 

cTq = 1 (70) 

are applied to Equation ( 68 ), the following results: 


J 

m 


(s) 




(71) 


or 


(72) 


If conditions of the standard uniaxial stress relaxation 
test, i.e.. 


a(t) = E^(t) 


where 




(73) 

(74) 


are applied to Equation (67), the following results: 


E„,(s) 


m 


or 






(75) 


(76) 


Combining Equations (72) and (76) yields 


E(n) = -5 

™ s^Jjjj(s) 


( 77 ) 
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since, from Equations (67) and (68), 




_1 


(78) 


From the definition of Poisson's ratio, v, as a function 
of time, 


e-T (t) 
v(t) = - 


(79) 


where e,j, is transverse extensional strain and is 
longitudinal extensional strain in a uniaxial test. If 
the uniaxial strain is held constant. Equation (79) becomes i 


£ 'T (t ) 

v(t) = - ^ 


(80) 


Taking the Laplace transform of both sides of this equation 
results in the following: 


_ ^ 'j> (s) 

v(s) (81) 

e 

o 

and therefore, 

(s) 

^ ■ I ( 82 ) 

(s) 


In the Laplace domain, Poisson's ratio is defined as: 

e-P (s) 
v(s) = 


( 83 ) 
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With the assximed constant uniaxial strain, 
Equation (83) becomes: 


v(s) * - 


e^(s) 


(84) 


v(s) = - 


t^^^(s) 


and therefore, 


^o 


SCip (s) 
v(s) 


Setting the right-hand sides if Equations (82) and (86) 
equal to each other, the following relationship is obtained: 


v(s) 


v(s) 


This equation reduces to 


v(s) 


_s 

v(s) 


v(s) = s v(s) 


^ich is the "associated elastic" Poisson's ratio as a 
function of the Laplace parameter. 


8.2 DETERMINATION OF CREEP COMPLIANCE, Jjj^(t) 

The creep compliance of a linearly viscoelastic material 
is determined by a curve fit to experimental data obtained 
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from a standard creep test. If a least squares curve fit 
is utilized, the creep compliance may be expressed 

in a Prony series formulation as [16]. 

m 

* D + V F exp (-R.t) + Ht (90) 

m i«l ^ ^ 


where D, and H are constants evaluated from the least 

squares curve fit to experimental data and the values 
are assigned. 

8.3 DETERMINATION OF RELAXATION MODULUS, E^^^Ct) 

The relaxation modulus and creep compliance are related 
by their Laplace transforms as given in Equation (78) . 

By inverting this relations.iip through the use of partial 
fractions, the expression for the relaxation modulus, E^(t) 
can be obtained. 

The Laplace transfoirm of the creep compliance in 
Equation (90) is: 


m 


J^(S) r 


F. 

1 


i=l s + R_. s 


(91) 


If m = 4, for example, the expanded form is 

F 


J (s) = - + 


Fi F2 F3 


m 


S+Rn 


s Rr 


s+R, 


4 . H 

+ 

s+R, 


which can be written as 


J_(s) 


C.s^ + C2S^ + + CjS + Cj 


m 


(92) 


(93) 


s^(s+R 3^) (sfR2) (s+R3)(s+R^) 

where the C^, i»l,2,...,6, are coefficients evaluated by 
taking a common denominator. 
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From the relationship given in Equation (77), 


m 


s\(e) 


(s+Rj^) (S+R 2 ) (s+R^) (s+R^) 

+ €28^ + €38^ + + C5S + Cg 


Solv5.ng for the roots of the fifth order polynominal in the 
denominator, and denoting the roots as Xj . where j » 
yields the expression: 


(s+R ) (s+R ) (S4K3) (s+R ) 

f„(s) 1 2 3 i_ 

(s+Xj^) (s+X2)(s+X3) (s+X^) (S+X3) 


(95) 


Inverting this equation by the method of partial fractions 
[ 17 ] gives the form 

e (s) - (96) 

® Q(s) 

where Q(s) and P(s) are polynomials and the degree of P(s) 
is less than the degree of Q(s). Equation (95) may be 
written as the sum of partial fractions: 



(s) 



s+X, s + X,, s + X« s + X, s+X 

12 3 4 


(97) 


where B ^ , j= 1,2 5, is evaluated: 

P (X.) 

B = J — (98) 

J Q' (Xj) 

The variable Q'(Xj) is the first time derivative of the 
polynomial Q (Xj) and P (X^) is the polynomial evaluated 
at the roots X^ . The resulting equation for the relaxation 
modulijs in the time domain is obtained: 
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5 

E„(t) » E B. exp (-A^t) (99) 



The bulk modulus for the epoxy resin material system is 
taken as a constant. This assumption is consistent with 
experimental evidence and accepted practice. The 
time -dependent expressions for the matrix properties, 
E^(t) and Vjjj(t) are thus determined. 


CHAPTER 9. ASSOCIATED ELASTIC SOLUTION 


In order to study the time -depen dent behavior of a 
filamentary composite spherical pressure vessel subjected 
to a constant internal pressure, the correspondence princi- 
ple was employed. Based upon this principle, the problem 
is transformed from the time domain to the Laplace domain, 
and the time -dependent problem is replaced by an "associated 
elastic" problem with constitutive relationships and bound- 
ary conditions expressed as ftinctions of the Laplace 
transform parameter s instead of the time parameter t. 
Inversion of the associated elastic solution then yields 
the solution for the time -dependent linearly viscoelastic 
problem. 

The basic inversion process can usually be accomplished 
by the smthod of partial fractions if the mathematical nK>del 
has a closed form elastic solution and simple material 
parameters [18]. However, for a typical actual material 
system, the function to be inverted is often known only 
for discrete positive real values of the transform parameter 
and, therefore, numerical techniques are appropriately used 
to obtain the transformed solution. Although nmaerical 
methods were not used to obtain the solution presented 
here, the complexity of the employed equations warranted 
the use of an appropriate numerical inversion technique. 
Aim>ng numerous techniques, the collocation method [19] is 
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apparently an accurate Inversion technique which yields 
a relatively straightforward calculation, regardless of 
the complexity of material representation. 

This approximate -inversion method yields a 
time -dependent sol'ttlon of the form [20]: 


Pf(t) - A + Bt + Z h^ e ®v (103) 

v=l 

where p£ is burst pressure and A, B, h^, and are 
constants. The constant A is evaluated from initial 
conditions, 

k 


A = P£ (t ) - I h 

^ ° vl ^ (104) 

In this case, an internal pressure p is applied at time t^ 
and held constant. Since the long-time value of p^ is 
assumed to be a finite constant, the value B * 0 substituted 
into Equation (103) yields: 


k t 

Pf(t) - A + E h^e' “v (105) 

V - 1 

Taking the Laplace transform of this equation and rearranr^ing 
gives : 


P£ (s) - A + s Z 


v*l I (1/ot^) + 


(106) 


where w ■ 1,2, ...,k. 
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Substituting Equation (104) Into Equation (106) and 
rearranging yields: 


k 

Z 

v-1 


(1 + -^ ) 


Pf (Cq) - 8 pf(3) 


(107) 


Because a six-term exponential function will be used to 
define the burst pressure p£ as a function of time, six 
discrete values of the Laplace parameter are utilized. An 
appropriate range of these values Is chosen. 

The selected expression for Is [20l: 


and 


■ e (v= 1,2,.. .,6) 

Equation (105) can be expressed as: 


(108) 


1 

1 

1 

2 

1+e"^ 

1+e’^ 

1 

1 

1 

l+e^ 

2 

1+e'^ 

1 

1 

1 

1+e^ 

1+e^ 

2 


1 

1 

1+e^ 

1+e^ 

1+e^ 

1 

1 

1 

IT? 

1+e^ 

1+e^ 

1 

1 

1 

l+el® 

l+e« 

1-^e^ 



T 


1+e 

1 


1+e 


1 

1 

1+e’® 

l+e’’^® 

1 

1 

1+e’® 

1+e’® 

1 

1 

1+e’^ 

l+e-^ 

1 

1 

1+e’^ 

l+e’^ 

1 

1 

2 

1+e"^ 

1 

1 


1+e^ 2 



Pf(to)-sPf(«)|,.e-5 

Pf(to)-8Pf(»)|,.e’^ 

Pf(to)-SPf(«)|g.e-l 

Pf(to)-8Pf(s) 

Pf(t^)-spj(s) lj.^3 

;Pf(to)-spj(s)|,.^5 


(109) 

Utilizing Equations (76) and (89) for the Laplace tranforms 
of the relaxation modulus and Poisson's ratio, respectively, 
expressions for s ^qj(s) and s v(s) are obtained. 
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These expressions are evaluated for discrete values of 
the Laplace parameter and inserted Into the constitutive 
equations for the laminate given in Chapter 6. The laminate 
properties In the Laplace domain are then used to calculate 
burst pressure using the maxlmiun strain theory of failure 
presented In Chapter 7. The associated elastic solution 
yielding p^Cs) Is thus obtained: 

Pf(B) - «K( , (110) 

+62n(s) 

f(8) 


The values of P£(s) and the burst pressure at 

time zero, are substituted into Equations (107) from which 
the h^ values arc computed. The value of A Is then calculated 
from Equation (104). Substituting these values of h^ and A 
into Equation (103) yields the expression for the time- 
dependent burst pressure for the filamentary composite 
spherical pressure vessel. 

A similar associated elastic solution may be used to 
calculate the critical strain, Eqq , as a function of time. 

The Laplace transform of a constant Internal pressure is 
used to calculate sT.. (s) : 


ee 

+ 32n(«) 




L. 

8 


f(8) 


g(s) 


(Ul) 
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The collocation n^thod of Inversion as outlined above Is 
then used to obtain the time -dependent strain at the 
Inner surface of the filamentary coai|)08lte spherical 
pressure vessel. 

The mathematical model for a linearly viscoelastic 
material therefore yields expressions vdilch describe the 
time -dependent response of a fllairantary composite struc- 
ture. A numerical example will be presented to show the 
application of this method to a Kevlar/epoxy spherical 
pressure vessel subjected to a constant Internal pressure. 


CHAPTER 10. APPLICATION 


The material system chosen to demonstrate the analysis 
technique presented is a Kevlar/HBRF>55 filamentary composite. 
The HBRF-55 matrix material is a resin currently used in 
the filas^nt wound motorcase of the Space Shuttle's Solid 
Rocket Boosters. The Kevlar fiber properties were supplied 
by the manufacturer and obtained from the literature [4] 
and are listed in Table 1. For the analysis, the fiber 
volume fractions (ratio of fiber volume to total laminate 
volun^) used were 0.5. 0.6, and 0.7. These reflect valueb 
attainable for the system when wound into a spherical geom- 
etry. To demonstrate the effect of increasing the compos ^ 
shell thickness, the solution was obtained for values of 8 
(ratio of outer radius to inner radius) of 1,1, 1.3, and 1.5. 

10.1 VISCOELASTIC MATERIAL PROPERTIES OF HBRF-55 RESIN 
As described in Chapter 8, a creep (or equivalent) test 
is necessary in order to determine the time -dependent 
properties of a linearly viscoelastic material . Using the 
apparatus shown in Flgute 4, a constant load was applied 
to a rectangular sample c* HBRF-55 resin. The deformation 
of the sample was c^asured by a dial indicator, and at 
selected time Intervals, the deformation of the sample 
was recorded. 


43 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


45 



eep tes 











46 


The strain was calculated from Equation (111) and 

these values are listed In Table 2. 

The FORTRAN program CREEP listed In Appendix A was written 

to calculate the creep compliance for each time Increment 

and to obtain a creep compliance expression as a function 

of time by applying a least squares curve fit. A 

2 

cross-sectional area of 0.258 In and an applied load of 
107.9 pounds, provided a constant uniaxial stress of 574.2 
psi. The creep compliance was, as outlined in Section 8.2, 
computed In program CREEP using Equation (90) for each 
time Increment. A least squares curve fit of these data 
yielded the following expression for Jjj^(t) , the creep 
compliance of the HBRF-55 epoxy resin matrix: 

Jj^(t) = 0.2979 X 10“^ + 0.2524 x 10"^ 

-0.22 X 10'^ e“^ -0.5191 x 10'^ e"^®^ 

-0.181 X 10"^ + 0.223 x 10"\ (112) 

A plot of the curve fit and the Input data points Is given 
In Figure 5. These plots versus log time are given in 
Figure 6. 

Following the development of equations in Sections 8.3 
and 8.4, the expressions for relaxation modulus and Poisson’s 
ratio as a function of time were obtained. The Laplace 
transform of is and is programmed into CREEP. 

The expression for Ejj^(s) is obtained from Equation (99) 
and a fifth order polynomial results in the denominator. 
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TABLE 2. EXPERIMENTAL CREEP TEST 

Sample Length «■ 1.952 inches 
Constant Applied Stress - 570.2 psi 


time 

LOG TIME 

AL 

GC 

(hrs . ) 

(hrs. ) 

(in. X 10"^) 

(in/ in x 10 

0.001 

-3.0 

31.0 

1.5881 

0.001778 

-2.75 

31.5 

1.6137 

0.003167 

-2.5 

32.0 

1.6393 

0.005611 

-2.25 

32.2 

1.6496 

0.01 

-2.0 

32.6 

1.6701 

0.017778 

-1.75 

32.9 

1.6855 

0.031623 

-1.5 

r 33.2 

1.7008 

0.05556 

-1.255 

33.5 

1.7162 

0.1 

-1.0 

33.8 

1.7316 

0.1778 

-0.75 

34.0 

1.7418 

0.31611 

-0.5 

34.5 

1.7674 

0.5622 

-0.25 

34.8 

1.7828 

1.0 

0.0 

35.2 

1.8033 

1.945 

0.289 

35.9 

1.8391 

3.1622 

0.5 

36.2 

1.8545 

5.79 

0.7627 

36.3 

1.8596 

10.0 

1.0 

37.0 

1.8955 

17.7828 

1.25 

38.3 

1.9621 


CURVE EIT 

EXPERIMENTAL DATA 



FIGURE 5. CREEP COMPLIANCE VS. TIME 
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COMPLIANCE VS. LOG TIME 



Roots of this polynomial were extracted using UNIVAC 
subroutine RTPOLY. The method of partial fractions Is 
utilized in CREEP to invert the Laplace expression 
Into the time domain. The relaxation modulus, Egj(t), Is 
thus obtained for the epoxy resin HBRF-55: 


E^(t) - 0.2212 X 10^ g-1065t ^ 0.5835 x 10^ g-lO.lSt 
+ 0.2242 X 10^ - 0.3064 x 10^ 

+ 0.3402 X 10® e'.0075t 

A plot of Ejj^(t) versus log time is given in Figure 7. 

The relaxation modulus of HBRF*55 at time zero was calcu- 
lated to be 360,000 psi. 

The expression for Poisson’s ratio, (t) » was obtained 
using the expression above for the relaxation modulus, E^(t) . 
Substituting into Equation (102) a bulk modulus for HBRF-55 
of 400,000 psi and rearranging, the equation for the Poisson' 
ratio is: 




1 

2 ■ 2,400,000 


(114) 


From the expression for E^(t) given in Equation (113) the 
resulting expression for v^jj(t) is: 

v^(t) » -9.2167 X 10"^e‘^°®*^*^ - 1.4313 x 10‘^e"^°* 
-9.3417 X lO’^e'^’®^^ + 1.2767 x 10"^e"*®^^*^ 
-1.4275 X 10“^e"-®°^^^ + 0.5 (115) 


400000.00 
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FIGURE 7. RELAXATION MODULUS VS. LOG TIME 
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A plot of v^(t) versus log time is shown in Figure 8. The 
Poisson's ratio for HBRF-55 at time zero was calculated to 
be 0.35. 


10.2 ASSOCIATED ELASTIC SOLUTION 

The correspondence principle, along with the collocation 
method, as derived in Chapter 9, may be used with the 
maximum strain theory of failure developed in Chapter 7 
to predict the burst pressure as a function of time for 
a filamentary composite spherical pressure vessel subjected 
to a constant internal pressure. A FORTRAN program, THESIS, 
was developed to generate the time -dependent response of 
a Kevlar/HBRF-55 system and is listed in Appendix B. The 
program predicts burst pressure and critical strain as a 
function of time. 


10.2.1 PROPERTIES OF THE LAMINATE 

The previously derived expressions for and 

Vm^t) are transformed into the Laplace domain utilizing 
Equations ( 76 ) and (89) and result in the equations: 


22,120 s 5,835 s 22,420 s 

sE (s) * + + 

“ s + 106.5 s + 10.18 s + 1.07 

- 30.640 s ^ 340,200 s 
s + .092 8 + .0075 

/.X „ c 0.0092 s 0.0024 s 0.0093 s 

sv_(s; * 0.5 - - - 

“ 8 + 106.5 s + 10.18 8 + 1.07 

+ 0.0128 s - 0.1418 8 


(116) 


(117) 


s + .092 


s + .0075 



FIGURE 8. POISSON'S RATIO VS. LOG TIME 
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The matrix properties and are thus calculated In the 
program THESIS and summarized In Table 3. These values are 
then substituted Into the lamina and laminate equations 
developed In Chapters 4, 5, and 6. For the case of 
V£ - 0.6 and 6 * 1.1, lamina properties computed In THESIS 
as a function of the six values of the Laplace parameter 
are given In Table 4. Laminate properties are listed In 
Table 5. Also listed are the computed values for n, the 
degree of anisotropy. For an Isotropic material system, 
n » 1.5, and the departure from this value Indicates the 
degree of anisotropy In the system. 

10.2.2 MAXIMUM STRAIN THEORY OF FAILURE 

The burst pressure and critical strain, utilizing 
Equations (61) and (59), respectively, are based upon the 
maximum strain theory of failure. In the Laplace domain, 
these expressions are given by Equations (110) and (ill). 

The ultimate clrcumferenctlal strain for Kevlar fiber Is 
0.0192 [ 4] and Is the numerical value for K^. 

Equation (105) In the form of Equation (109) Is then 
programmed into THESIS. The vector (PfC^Q) * sp£(s)} 
is obtained by calculating Pf (^ q) from Equation (61) in 
program THESIS. CCC, given in Appendix B. The value for 
sp£(s) Is subtracted from PfCt^) and evaluated for each 
discrete value of the Laplace parameter. From Equation 
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TABLE 3. VALUES OF a, E„(s), and v (s) 

01 m 


8 


''mCJ 

148.4 

350.2 

0.3542 

20.09 

338.2 

0.3591 

2.718 

327.5 

0.3636 

0.368 

314.9 

0.3688 

0.050 

286.1 

0.3808 

0.0067 

158.6 

0.4339 
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TABLE 4. VALUES FOR LAMINA PROPERTIES 


III 8 

ElC.) 

(v. “ 0.6, 

6 - 1.1) 

V li>(8) 

ijiip ( 8 ) 

5 

11.54 

0.8519 

0.5473 

0.2617 

0.3374 

3 

11.54 

0.8236 

0.5283 

0.2636 

0.3422 

1 

11.53 

0.7981 

0.5112 

0.2654 

0.3465 

-1 

11.53 

0.7682 

0.4913 

0.2675 

0.3517 

-3 

11.51 

0.6994 

0.4456 

0.2723 

0.3640 

-5 

11.46 

0.3916 

0.2453 

0.2936 

0.4251 
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TABLE 5 . VALUES FOR LAMINATE PROPERTIES AND n 
(v£ ■ 0.6, 3 ■ 1*1) 


In s 

"Cs), 

Crr<.) 



CerO 

Ej(s) 

Ej.{8) 

j 

5 

3.316 

0.9712 

5.118 

1.569 

0.3379 

4.579 

0.9371 


3 

3.356 

0.9424 

5.093 

1.569 

0.3319 

4.553 

0.9093 


1 

3.395 

0.9163 

5.071 

1.569 

0.3263 

4.529 

0.8842 

f 

-1 

3.442 

0.8855 

5.045 

1.569 

0.3194 

4.501 

0.8546 


-3 

3.564 

0.8141 

4.984 

1.567 

0.3026 

4.437 

0.7862 

1 

-5 

4.469 

0.4814 

4.697 

1.537 

0.2043 

4.154 

0.4680 

j 
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(109), the values for are calculated in double precision 
with UNIVAC subroutine DPINV. 

The constant A is evaluated from Equation (104), using the 
initial burst pressure and the sum of the values of h^. This 
operation is easily accoaq>lished in the program THESIS. 

Finally, the burst pressure in th-% tin» domain can be 
calculated from Equation (110), where the units cf time is 
hours. The burst pressures confuted from THESIS are plotted 
with respect to log time. Figure 9 shows burst pressure 
as a function of fiber volume fraction, Figure 10 

demonstrates the effect of varying g, the ratio of 
cooq>osite shell outer radius to inner radius. 

In a similar fashion, the vector {€QQ(t^)- 8 eQQ(s)}is 
obtained by calculating initial strain c (t^) in program 
THESIS. CCC. The value for se’ 00 (s) is subtracted from ^ggCt^) 
and evaluated for each discrete value of the Laplace 
paran^ter. Then, the values of h^ are calculated. The 
constant A' is evaluated and input into the time -dependent 
expression for strain. Computed values for strain from the 
program THESIS are plotted with respect to log time. 

Figure 11 Is a plot of strain versus log time and illus- 
trates the effect of fiber volume fraction, v^. Figure 12 
demonstrates the effect of varying g from 1.1 to 1.5. 
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A comparison of predicted burst pressure in Figure 9 
and strain In Figure 11 Is shown for fiber volume fractions 
0.5, 0.6, and 0.7. These plots show that an Increase In 
fiber volume fraction increases the performance of the 
composite. In particular, Figure 11 shows that the failure 
strain of Kevlar, 0.0192, Is approached almost immediately 
for V£ = 0.5, is approached much later for V£ » 0.6, and 
is never approached for V£ = 0.7. 

The effect of shell thickness on vessel performance Is 
shown by varying 8 from 1.1 to 1.5 in Figures 10 and 12. 

Figures 9 through 12 suimnarlze the results of THESIS, 
a program of the time -dependent burst pressure and strain 
equations from the associated elastic solution based upon 
the maximum strain theory. 


CHAPTER 11. CONCLUSIONS 


The maximum strain theory of failure developed for 
composite filamentary pressure vessels, when combined with 
the linearly viscoelastic behavior of the composite matrix, 
yields a solution which predicts burst pressure and critical 
strain for the vessels as a function of time. This solution 
is obtained by solving the equations in the Laplace domain 
and inverting them back into the time domain using the 
method of collocation. An evaluation of the burst pressures 
and strains computed using this method demonstrates that 
when the time-dependent response of a composite matrix is 
considered, the failure performance of the vessel is reduced. 
It is thus concluded that this general method may be used to 
predict a time -dependent response, and that the linearlv 
viscoelastic properties of a composite matrix contribute to 
this reduction in vessel performance. Knowing the relaxation 

modulus and Poisson's ratio of a matrix with respect to time, 
design criteria for a filamentary composite spherical 
pressure vessel could be derived which consider the 
time -dependent response of the matrix. 

The numerical evaluation of the Kevlar-HBRF-55 system 
illustrates the reduction in burst pressure and Increase in 
critical strain with respect to time. The inclusion of the 
viscoelastic material properties of the resin into the 
associated elastic solution is thus a viable approach for 
predicting time -dependent vessel performance. 
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There are additional mechanisms which affect vessel 
performance in time, such as nonlinear effects and/or 
matrix crazing, perhaps Initiated by this modeled time- 
dependent matrix stiffness reduction. The same type 
analysis could include a damage criterion or function and 
provide a more accurate prediction of time-dependent 
response. It has been shown that this type of analysis 
can be used to predict the time -dependent response of 
filamentary composite spherical pressure vessels and 
can possibly be used to identify the mechanism(s) which 
cause the observed failures. 
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